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ABSTRACT 

Binaries consisting of a pulsar and a black hole (BH) are a holy grail of astrophysics, 
both for their significance for stellar evolution and for their potential application as 
probes of strong gravity. In spite of extensive surveys of our Galaxy and its system 
of globular clusters, no pulsar-black hole (PSR-BH) binary has been found to date. 
Clues as to where such systems might exist are therefore highly desirable. We show 
that if the central parsec around Sgr A* harbors a cluster of ^ 25, 000 stellar BHs 
(as predicted by mass segregation arguments) and if it is also rich in recycled pulsar 
binaries (by analogy with globular clusters), then 3-body exchange interactions should 
produce PSR-BHs in the Galactic center. Simple estimates of the formation rate and 
survival time of these binaries suggest that a few PSR-BHs should be present in the 
central parsec today. The proposed formation mechanism makes unique predictions 
for the PSR-BH properties: 1) the binary would reside within ^ 1 pc of Sgr A*; 2) the 
pulsar would be recycled, with a period of ~ 1 to a few tens of milliseconds, and a low 
magnetic field B < 10^° G; 3) the binary would have high eccentricity, e ^ 0.8, but 
with a large scatter; and 4) the binary would be relatively wide, with semi-major axis 
Ob ~ 0.1— > 3 AU. The potential discovery of a PSR-BH binary therefore provides a 
strong motivation for deep, high-frequency radio searches for recycled pulsars toward 
the Galactic center. 
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1 INTRODUCTION 

Pulsars in binary systems provide accurate clocks that can 
be used to infer the properties of the binary orbit and its stel- 
lar components with precision. In addition to yielding im- 
portant constraints on stellar evolution, binary pulsars have 
been used to test the validity of general relativity and to put 
stringent constraints on alternative theories of gravity (see, 
e.g., the review article by Stairs||2004 1. The original binary 
pulsar PSR B1913+16 ( |Hulse fc Taylor|1975[ ), consisting of 
a pulsar in orbit around another unseen neutron star, per- 
mitted the first tests of general relativity in the strong field 
regime and the most compelling (albeit indirect) evidence 
for the gravitational wave emission, via the measurement 
of its orbital period derivative. Today, the best strong-field 
tests of general relativity are provided by the "double pul- 
sar" PSR J0737-3039 | |Burgay et al][2003l ), the onl y known 
example of two pulsars in orbit around each other ( Kramer 
fc Stairs|2008l ) 



fraction of pulsars should have black hole (BH) companions 



(e.g.,|Narayan et al.j l991[|Phinney|1991|[Po rte gies Zwart fc 
Yungclson 1998; Bcthe & Brown 1999, Sipior fc Sigurdsson 
2002^; Bclczynski ct al.,,2002, ,Pfahl et al.,,2005) and pulsar 



While theoretical considerations suggest that a certain 



black hole (PSR-BH) binaries are a holy grail of pulsar as- 
tronomy, none has been found to date. The discovery of such 
a system would represent a major step forward both from 
the astrophysical point of view and for the unique gravity 
tests that it might allow. Clues as to where PSR-BH sys- 
tems might be found are therefore highly desirable. Most 
theoretical studies have focused on PSR-BHs formed from 
primordial binaries in the Galactic field. As dense systems 
conducive to interactions between stars, globular clusters 
have also been proposed as sites where PSR-BH binaries 
may reside ( Sigurdsson|2003"l ) . Indeed, the high rates of stel- 
lar encounters in globular clusters explain the overdensity 
of recycled millisecond pulsar^ (MSPs) and X-ray binaries 
in globular clusters, by orders of magnitude, relative to the 
field (e.g., |Verbunt|1987[[Hu"t et al.|p9T] |Rasio et al!][2000 
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^ A recycled pulsar is one that has been spun up after birth by 
accretion of material from a binary companion. 
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Pooley et al.||2003| [Pooley fc Hutl|2006[ ). It is natural to ex- 
pect that analogous processes might lead to the formation of 
PSR-BH systems. An important limitation of this argument, 
however, is that BHs may be excessively rare in globular 
clusters. To date, there is no consensus on whether globular 
clusters host massive central BHs similar to those commonly 



found at the centers of galaxies (e.g., Gebhardt et al.|2005 
Noyola et al.„2008| [van der Marel fc Anderson||2010[ ). Fur- 



thermore, stellar-mass BHs are also expected to be few in 
globular clusters as a result of dynamical interactions that 



eject them (Sigurdsson & Hernquist 1993 Kulkarni et al 



1993[ [Kalogera et al.|2004| [Moody fc Sigurdsson|2009| ). 

In this paper, we consider the formation of PSR-BH 
binaries in another kind of dense stellar environment, the 
center of our Milky Way. Indeed, the nuclear star cluster at 
the center of our Galaxy should allow many of the stellar 
interaction processes taking place in globular clusters to op- 
erate, and in addition may be particularly rich in BHs. In 
fact, the surrounding bulge provides a large reservoir of stel- 
lar remnants. Because stellar (~ 10 M©) black holes are the 
most massive components of old stellar populations, they 
sink to the center owing to dy namical friction (|Morris|1993[ 



Miralda-Escu de fc Gou ld|2000|[FYeitag et al.|2006| ). |Miralda; 
Escude & Gould (2000) showed that about 25,000 steUar are 



expected to have migrated into the central parsec in this 
way. Here, we show that these BHs are likely to form bi- 
naries with pulsars via exchange interactions, and that sev- 
eral MSP-BH binaries should survive in the central parsec 
today. This makes the Galactic center (GC) a promising ex- 
perimental target not only for pulsars orbiting the central 
supermassive black hole (Sgr A*;[Paczynski fc Trimble|1979[ 



Cordes fc Lazio|1997| [Pfahl fc Loeb[[2004[ ), but also stellar- 
mass BHs. Similar processes should operate in the nuclei of 
other galaxies, and may actually be more efficient in low- 
mass galaxies with no or low-mass central black holes (e.g.. 
Miller fc Lauburg|2009 1. Because of its proximity, the center 
of the Milky Way is however the most interesting observa- 
tionally at present. 

No pulsar has been found in the central parsec so far 
| |Deneva et al.|2009[ ). This dearth of observed pulsars in the 
GC is understood to arise from the strong scattering of radio 
waves by the dense, turbulent, ionized plasma in that region 
(e.g., [Lazio fc Cordes[[1998| [Johnston et al.[ [2006'), rather 
than by an intrinsic absence of pulsars. If the scattering of 
radio pulses broadens them by a time scale comparable to 
or longer than the pulse period, the signal becomes effec- 
tively constant and the pulsar disappears. Since the scat- 
tering broadening time scales with frequency as Tscat oc 
(e.g., [Cordes fc Lazio[[1997[ ), radio pulsars can nevertheless 
in principle be detected in the GC by observing at suffi- 



ciently high frequency. Recently, Macquart et al. (20101 re- 



ported on a 15 GHz search for radio pulsars in the central 
parsec with the Green Bank Telescope (GBT), the highest 
frequency search for pulsars toward the GC to date. Obser- 
vations at this frequency are sensitive to pulsars with peri- 
ods > 50 ms. Although it was not confirmed in subsequent 



data taken in 2008, Macquart et al.[ ( [2010) 1 did detect a 607 
ms pulsar candidate in their 2006 observations. It must be 
noted that the 2008 non-detection does not rule out that the 
candidate is a real pulsar, as a source in a short-period orbit 
(<100 yr) orbit around the GC could have its emission beam 
rotated away from our line of sight on a 2-yr time scale. Al- 



ternatively, the pulsar could have moved sufficiently in the 
inhomogeneous GC plasma to induce significant variations 
in its pulsed flux, owing either to scattering or scintillation. 
Regardless of the reality of this pulsar candidate, the upper 
limit of 90 ordinary pulsars in the central parsec inferred by 



Macquart et al. ( 2010 1 is comparable to the expected num- 



ber based on the simple expectation that there should be 
about as many ordinary pulsars as progenitor massive stars 
of mass > 8 M0 (e.g., Bartko et al.[[2010 l, since their life- 
times are comparable. Further searches are therefore highly 
warranted, as even current technology stands a good chance 
of detecting radio pulsars in the central parsec. 

For our purpose, it is especially noteworthy that even 



the Macquart et al. (20101 15 GHz periodicity search, and 



by extension all previous searches at lower frequencies, were 
completely insensitive to MSPs with pulse periods P ~ 
1— few ms. Existing observations thus would not have de- 
tected the MSP-BH binaries that we predict. Nevertheless, 
pulsars have been detected at a variety of higher frequen- 



cies, from 32 GHz (e.g., Lohmer et al. 20081, to 43 GHz 



(e.g., [Kramer et aLl[1997[), to 87 GHz^(e.g ., [Morris et al. 
[Tgffl ), to 144 GHz (e.g., [Camilo et al.[[2007[ ). Furthermore, 
the radio spectra of MSPs are similar to those of ordinary 
pulsars (e.g., Kramer et al.|1998 Toscano et al.[[l998 l, sug- 
gesting that radio observations at frequencies high enough 
to beat down interstellar scattering in the GC for millisec- 
ond periods are possible. In theory, GC MSPs could also be 
observed in higher-energy bands, including in the X— rays 
and 7— rays. 

In the rest of this paper, we describe our new formation 
scenario for MSP-BH binaries in the Galactic center and 
present simple estimates of the formation rate and survival 
of these systems ([j2|. Because of the technical challenges of 
simulating the dynamics of the GC, including the critical bi- 
nary processes and the effects of different stellar masses (e.g., 
Freitag et al.|20d6 Hopman|2009 1 , and since substantial un- 
certainties exist, we focus here on basic analytic considera- 
tions. Our here is to outline the relevant physical processes 
and to combine them in order-of-magnitude estimates for 
the number of MSP-BH binaries that should survive in the 
GC today, as a motivation for observational efforts as well 
as further theoretical studies. In ij3| we discuss the prospects 
for detecting the predicted systems, the unique signatures of 
the proposed formation channel, and the potential implica- 
tions for physics and astrophysics. 



2 MSP-BH BINARIES IN THE GC 

The MSP-BH binary formation scenario that we envision 
in the GC is as follows. We assume that, in analogy to the 
phenomenon in globular clusters, the high stellar densities 
in the central parsec lead to the dynamical formation of 
MSPs. We further assume, again by analogy with globular 
clusters, that most of these MSPs will be found in binaries 



with white dwarf (WD) companions. In [ 2.4 we discuss how 



the properties of globular clusters are scaled to the nuclear 
star cluster. What is unique about the GC is that the MSP- 
WD binaries have a significant probability of undergoing a 
3-body interaction with a stellar BH in the central cluster. 
A likely outcome of many such interactions is the exchange 
of the WD companion for the stellar BH, yielding a MSP- 
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BH binary. In what follows, we investigate this scenario by 
quantifying the time scales for the formation and survival of 
the MSP-BHs. 

To simplify the calculations, we assume a steady-state 
stellar background and calculate the rates of MSP-BH for- 
mation and destruction in this background. This approxima- 
tion is justified by the typical lifetimes less than a few Gyr of 
the MSP-BH binaries (^2.3 and ^2.5l, i.e. significantly less 



than the age of the Galaxy and therefore presumably of the 
nuclear cluster. Potential effects arising from recent events 
are outlined in 93.41 



2.1 Stellar Distribution in the GC 

The supermassive black hole at the center of our Galaxy, 
Sgr A*, is taken to be at the origin and to have a mass 



Msi. 



4x10" Mq (e.| 



Ghez et al. 



2008 



Gillessen et al. 



20091. Schodel et al. (20091 examined the proper motions of 
stars out to a distance of 1 pc from Sgr A* and estimated the 
extended mass (excluding that of the SMBH) to be 0.5— 1.5x 
10® M© within that radius. We assume that this extended 
mass follows a power-law cusp 



ro 



(1) 



normalized such that M*(< Ipc) — 4tt J^^'^ drr p*{r) = 10 
M0 . The power-law index 7 is not uniquely constrained ob- 



servationally ( Schodel et aL||2009 l but consistent with the 
value 7 = 1.3 predicted by IFreitag et al. (20061 for main 



sequence stars in their standard Milky Way nucleus simula- 
tion. We adopt this value in our calculations, but note that 
it has a small impact on our results since the SMBH dom- 
inates the velocity dispersion throughout most the region 
of interest (r < 0.5 pc). The gravitational potential of the 
combination of the SMBH and the extended central mass 
distribution is then 



<l>(r) - <l>(ro) 



where 



tiro) 



tir) GMs^ 



2 I 

Wccxt I 



r) 



AixGp.. 



.o'"o , 



.2-7 



(7/2), 

(2) 
(3) 



(4) 



For a spherical system with an isotropic velocity distribution 
at each point, the Jeans equation reads 

d(ncr2) _ _ d<E> 
dr dr 

where n is the number density of the constituents of interest 
and a is their velocity dispersion (e.g., [Blnney fc Tremaine 
20081. For the potential in equation ([2| and a stellar popu- 
lation following n oc r~", this yields 

47rGp*,or|]' 



GAfsMBH 

(1 + a)r 



(5) 



(3-7)(7 + a-2) 

If gravitational encounters between stars are sufficiently 
frequent, they efficiently exchange energy. In equipartition, 
this would imply that more massive stars should have a 
lower velocity dispersion than less massive ones by a fac- 
tor ^JM\JM2, where the Mi are the masses of two stellar 
species under consideration. Self-gravitating systems, how- 
ever, are subject to the counter-acting effect of mass segre- 
gation, which tends to bring the more massive components 



inward, thereby leaving them on higher-velocity orbits (e.g., 
'Khalisi et al. 20071. Within the radius of influence of the 



SMBH, the actual squared velocity dispersion of a given 
species scales instead with the power-law index of its num- 
ber density distribution (the [1 -I- a]~^ factor in the second 
term of eq. ([5|). The dispersion of the relative velocities be- 
tween any two species, relevant for interaction cross sections, 

is then simply iT*,rei ~ \J'^'i 1 + 2 • 

Because the relaxation time is short at the Galactic 
center, old stars are expected to have settled in a cusp of 



approximately power-law form (e.g., Peebles 1972 Bahcall 



& Wolf 1976 19771. We therefore use simple power laws 



to describe the spatial distributions of stars. Of particular 
importance for our problem, stellar BHs of mass ~ 10 M0 
from the surrounding bulge should sink into the central par- 
sec on a dynamical friction time scale, which is about 10 
Gyr (roughl y the age of the Galaxy) at a radius of 5 p c 



from Sgr A* (Morris 



1993 



Miralda-Escude & Gould 2000k 



This process acts to concentrate ~ 25, 000 stellar BHs in 
the central parsec. Since it takes « 30 Gyr for these BHs 
to be swallowed by Sgr A* ( |Miralda-Escud6 fc Gould|2000[ ), 
most of these stellar BHs should remain there today. We 
parameterize the number density distribution of stellar BHs 
by 



nBH(r) = riBH.o 



(6) 



where we conventionally deflne ro = 0.5 pc (characteristic 
of the size of the stellar BH cluster) and the constant nen.o 
is set by normalizing to the total number of BHs, A^bh- The 
other stellar population of greatest interest for us are MSPs, 
whose distribution we parameterize as 



nMSp(»') = nMSP.o 



(7) 



For a single-mass stellar population, [Bahcall fc Wolf] 



( 1976 1 showed that the cusp approaches a universal slope 
of —7/4. In a more realistic stellar population consisting 



of multiple masses, Bahcall & Wolf (19771 argued that the 



slope would depend weakly on the mass of the component 
of interest, a — 3/2 + p, where p « 0.3Mi/Mi and Afi is 
the mass of the heaviest stellar component. In our fiducial 
calculations, we assume that the MSP- WD spatial distribu- 
tion follows that of neutron stars, and set qbh ~ 1-75 and 
Qmsp ~ 1-3, also consistent with the standard Milky Way 



nucleus simulation of Freitag et al. ( 2006 1 



Two variations about these fiducial choices warrant spe- 
cial consideration. First and most importantly, some obser- 
vations suggest that the stellar distribution near Sgr A* may 
not be cuspy but instead form a core or even a cavity within 
~ 0.2 pc ( |Do et al. 12009] [Buchholz et al.|2009l ). These obser- 
vations however currently only directly probe late-type red 
giants, and so do not necessarily exclude the presence of a 
relaxed cusp of stellar remnants as predicted theoretically]^ 



^ As the stellar BHs migrate inward, they eventually dominate 
the stellar distribution and their further evolution is determined 
by interactions with each other, rather than dynamical friction 
(e.g., [Freitag et al.|2006^ . 

For instance, physical collisions could preferentially deplete gi- 
ants from the inner cusp (e.g., Dale et al.|2009| ). 



© 0000 RAS, MNRAS 000, 000-000 



4 Claude- Andre Faucher-Giguere and Abraham Loeb 



Nevertheless, we will explore this possibility by truncating 
the PSR-BH formation rate at minimum radii i?min = 0.1 



and 0.2 pc (S 2.5 I, and show that our results are only weakly 
dependent on the presence of an inner cusp. Second, there 
is the "strong" segregation regime discussed by lAlexander] 
& Hopman (20091, in which the most massive components 



contribute only a minor fraction of the gravitational poten- 
tial and concentrate more strongly owing to dynamical fric- 
tion off less massive components in the stellar cusp (see also 



Keshet et al.|2009||Preto fc Amaro-Seoane|2010||Preto|2010[ 



Amaro-Seoane &: Preto||2010 K It is possible that the Milky 
Way nucleus marginally satisfies the criteria for strong seg- 



regation ( Hopman fc Alexander, 2006| ) , in which case a value 
as steep as obh ~ 2 could be realized. Our simple arguments 
are unfortunately not well suited to treat this singular case, 
in which the number of BHs inside a radius r is divergent for 
«BH 5S 2 if we do not account for an inner turn over. Further- 
more, recent results suggest that resonant relaxation could 
carve out a cavity similar to that observed in late-type stars 
in the inner ~ 0.2 pc ( [Madigan et al. 20101, so that the 
presence of such a singular inner cusp is not at present well 
supported for the Milky Way. We will therefore not consider 
this possibility further in this work. 



2.2 MSP-BH Formation Rate 

Most of the theoretical work on 3-body encounters has fo- 
cused either on the very soft regime (in which the binding 
energy of the binaries is much smaller than the kinetic en- 



ergy of background stars; e.g.. Hut 1983b I, the very hard 



regime (in which the binding energy of the binaries greatly 
exceeds the kinetic energy of the other stars; e.g. |Heggie| 



et al.|[T996 l, or on situations in which all the stars of have 
identical mass (e.g., |Hut fc BalicaIi|[T983l [Hut|[T984| . For 

our problem, many binaries straddle the soft-hard transi- 
tion and the principal stellar components (WDs, PSRs, and 
BHs) have substantially multiple masses, which has an im- 



portant impact on the reaction rates (e.g., Sigurdsson & 



Phinney 1993 1. To ensure accuracy, we rely on numerically 



evaluated 3-body cross sections that are valid regardless of 
the hardness of the binaries, and which incorporate the ef- 
fects of different masses. To do so, we use the sigmaS pro- 
gram from the publicly- available Starlab package ( |McMillan| 
fc Hut|1996 1. This program automatically performs a series 
of direct 3-body numerical experiments and determines the 
cross sections using a Monte Carlo technique. We then eval- 
uate the interaction rates (av) (where a denotes the cross 
section and v denotes the relative velocity between the bi- 
nary center of mass and the third star) by averaging over a 
Maxwellian velocity distribution at each radius from Sgr A* . 
The accuracy parameters of the automated sigma3 program 
and of the integration over the Maxwellian distribution are 
increased until convergence is attained, with negligible re- 
sultant statistical errors. For simplicity, we assume that all 
the WDs have a fixed mass Mwd, that all the neutron stars 
have a mass Mns, and that all the stellar BHs have a mass 

AfBH. 

For reference. Figure [l] shows 3-body interaction time 
scales. 



^cxch,ion (^) 



nBH(r-)(crcxch,ionf)(»") 



(8) 



for relevant combinations of stellar component masses. In 
both panels, the interactions are assumed to take place be- 
tween binaries and external stellar BHs of mass Mbh = 10 
Mq, with a power-law density profile of index cibh = 1.75 
and riBH.o = 10* pc"'' at rg = 0.5 pc, and for the nuclear 
mass distribution model described in the previous section. If 
the BH cluster were more concentrated in reality, the time 
scales would be correspondingly shorter, and vice versa. 

We first consider the spatially-averaged rate at which 
an individual MSP- WD system undergoes an exchange in- 
teraction yielding a MSP-BH system, for a fixed semi-major 
axis of the original binary, Ob (the "creation rate"): 

Cmsp-bh(A/wd, Mns, Mbh, a-b) = 



X 4tt 



drr nMSp{r)nBii{r) 



hv){MwD, Mns, Mbh, ab, 
< R 



(9) 



where AfMSp(-Rmin < r < i?niax) = 47r J^^^^^ drr^nmspir) 
is the total number of MSP- WD systems, and crj^^ is the 
cross section for an exchange interaction that ejects the WD 
and leaves behind a MSP-BH binary. We implicitly assume 
that all MSPs have a WD companion. Rmin and -Rmax corre- 
spond to the minimum and maximum radii at which MSP- 
BH formation proceeds. Very close to the SMBH, binaries 
are tidally separated, so that i?min > 0; we will however 
show that our results are weakly dependent on the exact 
value of -Rmin and fiducially adopt -Rmin ~ 0.1 pc. iimax sim- 
ply corresponds to the outer radius of the stellar BH cluster 
in which the exchange interactions take place. In their simple 



analytic modeling, ^Miralda-Escude fc Gould ( 2000 1 predict 
fimax = 0.7 pc, while the numerical calculations of |Freitag| 



et al. ( 2006 1 suggest that the BH density profile has already 
steepened substantially at r « 0.5 pc. We adopt this latter 
value for -Rmax in our fiducial estimates. 

As Figure [l] shows, the 3-body interaction cross sec- 
tions are sensitive to the binary semi-major axis. Because 
semi-major axes typically have wide distributions, it is im- 
portant to average the creation rate over such a distribution. 
We assume that the semi-major axis distribution is flat in 
logarithmic space between specified limits: 



dP 
lln Qb 



< flb < fflb ,1 







otherwise 



(10) 

This simple analytic form is a crude approximation to the 
observed distribution of binary pulsar semi-major axes and 
is a reasonable first-order approximation for our exploratory 
estimates at the GC, where there are currently no empir- 
ical constraints. We adopt (ab.min, ab,max) = (0.01, 0.3) 
AU in our fiducial calculations. The lower limit corresponds 
roughly to the orbital separation at which circular MSP- WD 
binaries merge in < 1 Gyr owing to the emission of gravita- 



tional waves (see eq. ( 19 1 below). The upper limit is roughly 
the maximum separation of MSP- WD binaries that survive 
evaporation for > 1 Gyr (Fig. [ijand S2.3I. 
We define: 



d In Qb 



dP 



■CMSP-BH(ab). (11) 



dlnflb 

In the absence of destruction processes and assuming a 
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Figure 1. Time scales for 3-body interactions of (single) binaries with a nuclear cluster of A/bh = ™3 = 10 Mq BHs, as a function of 
radius from Sgr A* and for different binary semi-major axes. Left: MSP-WD binaries with A/wD = "il = 1 Mq and AfjMS = ™2 = 1-4 
Mq. Right: MSP-BH binaries with MgH = = 10 Mq and Af^g = m2 = 1.4 Mq. These plots assume the stellar distribution model 
for the Galactic center described in i j2.l[ i.e. a central SMBH of mass AfgMBH = 4 X 10® Mq surrounded by a cluster with total mass 
M^,{< 1 pc) = 10® Mq within 1 pc and a power-law density slope 7 = 1.3. The stellar BHs are assumed to have a number density 
"BH = 10^ pc~^ at ro = 0.5 pc, and radial profile of slope obh = 1-75. For each process (line style), the time scale increases with 
decreasing semi-major axis. 



steady creation rate, the total number of MSP-BH systems 
formed after a time equal to the age of the Galaxy, tcai, 
would be: 



(C'MSP-BH)abiGalA^MSP- 



(12) 



For our fiducial parameter values and A^bh = 25, 000 stellar 
BHs between Rmin = 0.1 pc and i?max = 0.5 pc, 

-^Ms'p-BH ~; Q Q4 {CMSP-BH)ab ^ ^ ^Gal j ^^^^ 



Nmsp 



V{CMSP-BH)S^y \ W Gyr 



suggesting that up to a few percent of MSPs in the GC could 
be in BH binaries. 



MSP-BH binaries. Within the cluster of stellar BHs, binary 
ionization is dominated by encounters with the BHs owing 
to the ml dependence of the ionization cross section. 



where rris is the mass of the third star (e.g.. Hut 1983b I 



This is because the ionization rate scales as pams (where 
p3 = n^ma and ns is the number density), and this factor is 
largest by a factor of several for BHs in the standard Milky 



Way model of Freitag et al. ( 2006 1 on which we base our 
calculations. 

Ionization proceeds in two regimes. First, the binaries 
can be ionized in single encounters. In the limit of a high- 
velocity encounter, the time scale for this process is 



2.3 Survival Time of the MSP-BH Binaries 

In practice, the MSP-BH binaries are continuously de- 
stroyed. A more realistic estimate of the number remaining 
at the present time is therefore 



-^MSP-BH = {C'MSP-BHisurv)ab-^MSP, 



(14) 



where tsurv is the time for which a given binary survives. 



Equation (141 is easily understood as the number of MSP- 
BHs created less than a time tsuw in the past, and the sur- 
vival time can be expressed as the minimum time of the 
destruction time scales for the different relevant processes. 

We now summarize the different processes that can 
destroy the MSP-BH binaries in the GC and their corre- 
sponding time scales. 



80 



1/2 



0-3 



mi + m2 



ma 



(15) 



where we have defined tion = (?i3{o"ion«))~^ and 0-3 is the 
relative velocity dispersion between the third bodies and the 
centers of mass of the mi, m2 binaries. The cross section (Tion 



for ionization was derived by Hut ( 1983b I, and we averaged 



over a Maxwellian velocity distribution. We use a primed 
a(, to denote the semi-major axis of the MSP-BH formed in 
the exchange interaction, which in general differs from the 
semi-major axis of the original binary. 

The cumulative effect of weak encounters can also grad- 
ually evaporate the binaries. Generalizing the Fokker-Planck 



Ionization: Encounters with field stars can ionize the derivation in Binney & Tremaine ( 2008| ) for the case of equal 
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masses, the time scale for the evaporation process is 



mi, m2, and semi-major axis ab gets within a distance 



^cvaD — 



16 



1/2 



0-3 



77117722 



Ga'^pz in A 



777S 



(7711 + 7n3)(7n2 + TTia) 
277117712 + 7713(7711 + 77I2) 



where 



A; 



1/2 



2G[min(77ii, 7712)+ 7773] 



(16) 
(17) 

(18) 



For mi = A/bh, m2 = Mns, and ms = Mbh, appropriate for 
MSP-BHs in a cluster of stellar BHs, tevap ~ 0.4fion/lnA. 
For characteristic values a(, = 1 AU and (J12 = oz = 200 
km s~^, A ~ 7.2 and so fcvap ~ 0.2fion, only logarithmically 
dependent on the binary semi-major axis and surrounding 
velocity dispersion. The diffusive evaporation process 
therefore dominates over discrete ionization, and so we 
use a tidal separation time scale of fovap = 0.2tion in our 
calculations. This is a conservative approximation, as tevap 
is actually within a factor ~ 2 of Uo-n for the binaries that 
are only marg inally soft ( |Hut|1983a| |Rasio et al.|1995[ ). 



Exchange: When MSP-BHs encounter other stars, an 
exchange interaction may take place that ejects one of the 
binary components. As for ionization, this process is domi- 
nated by encounters with other stellar BHs, and in this case 
only the exchanges that eject the MSP effectively destroy 
the MSP-BH. The time scale for this process is evaluated 
as in equation ([8|, using cross sections numerically com- 
puted with Starlab. In most cases, exchanges are negligible 
relative to binary evaporation (see the right panel of Fig.jl]). 

Gravitational wave-dnven merger: A MSP-BH binary 
will merge owing to the emission of gravitational radiation 
after a time 



tGW « 200 Gyr 



0.1 AU 

/ Mns Mbh 
V 14 M| 



/ Mns + Mbh \ ^ 
V 11.4 M© J 



(1 



(19) 



where e' is its eccentricity of the MSP-BH formed 



in the exchange interaction (Peters 1964 Pfahl et al. 
2005). In the hard limit (where the exchange cross 
section is important), the semi-major axis of the new 
binary can be estimated as aj, ~ ab(MBH/MwD) 
( Sigurdsson & Phinney 1993 1 and the new binary ec- 
centricity is high, since 1 — e' ~ Rp/a'^, where Rp ~ 
ab(^wD/MBH)[(M wD + Mns + Mbh)/(Mwd + Mtis)f'^ 
( Heggie et al.|[l996l ). We use these analytic approximations 
for the semi-major axis and eccentricity in our calculations. 
In the limit Mwd ^ Mns ^ Mbh, the eccentricity ex- 
pression simplifies to 1 - e' ~ (Mwd /Mbh) (Mbh /Mns )^''^, 
indicating that the key parameter is Mwd /Mbh. For our 
fiducial masses, aj, ~ lOob and e' ~ 0.8. Because of the 
strong scaling of few with and since the semi-major axis 
is increased by a factor Mbh/Mwd ~ 10 in the exchange 
interaction, gravitational wave-driven merging is ultimately 
unimportant. 

Tidal separation: If a binary with component masses 



rr 



ah 



Ms 



1/3 



mi + m2 
f mi +m2 



3.3 X 10^^ pc 



( ) 
VO.l AU/ 



-1/3 



1/3 



(20) 



MsMBH 

Vll.4M0"y' ^4 X 10^ Mq^ 

of the central massive black hole, the tidal force will separate 
it. Tidal separation in the central parsec occurs predomi- 
nantly via the loss cone, i.e. through angular momentum 
diffusion which can shrink the periapsis distance from the 



SMBH by increasing the orbital eccentricity (e.g., Frank & 



Rees|1976|[Lightman fc Shapiro|197"7l|Yu fc Tremaine|2003 1 

A simple estimate of the time scale for separating MSP- 
BH binaries via the loss cone is obtained by noting that the 
process is analogous to the one responsible for the direct 
swallowing of stellar BHs by the SMBH, but with a larger 
destruction inner radius Qmin equal to rT. |Miralda-Escude fc| 
Gould (20001 showed that the time scale for depleting the 
stellar BH cluster by direct swallowing is ti^ ~ 30 Gyr and 
scales only logarithmically with qmin as 

In {Rma^/q mill ) 

In [(A/sMBH/i\fBH)(gmin/i?max)l/4] ' 

where i?max is the maximum integration radius, taken to 
correspond to the radius of the BH cluster. The factor in 



tic oc 



(21) 



equation (21 1 is reduced by less than 40% if we consider 
the tidal separation of MSP-BH binaries with Ob = 0.1 AU 
instead of direct BH swallowing, for iimax ~ 0.5 pc. The 
physical reason for the weak dependence on Qmin is that the 
process occurs in the "empty" loss cone regime, in which an 
orbiting mass is only deflected by a small amount during 
an orbital period. This allows it to effectively sweep over all 
possible eccentricities, until the narrow range corresponding 
to q^nin is entered. 

MSP-BH binaries are therefore not tidally separated 
in the loss cone for nearly as long as it takes to deplete 
stellar BH cluster. Since the time scale exceeds the age of 
the Universe, we neglect this process. 

Radial wandering: Our discussion has so far implicitly 
assumed that binaries evolve at fixed radii from the central 
SMBH. In reality, most binaries have eccentric orbits 
whose orbital elements change with time as they undergo 
encounters with other stars, which causes them to wander in 
radius. In our picture of a steady-state stellar background, 
this does not affect the average MSP-BH creation rate. 

The destruction time scale for any given system must 
however be weighted by the amount of time that the sys- 
tem spends at different radii. In the limit of a short wan- 
dering time scale, the cusp is well mixed, with a given 
MSP-BH spending a time in a shell of thickness Ar pro- 
portional to the average number of systems in this shell, 
47rr'^7iMSP-BH(»') Ar. This limit is reasonable, because the 
orbital time is negligible, and the two-body relaxation time 
scale for 10 Mq BHs is < 10® yr. We therefore expect 
the MSP-BHs to generally have time to explore a range of 
semi-major axes and eccentricities in their orbits around the 
SMBH. 

Given an ionization/exchange survival time for a MSP- 
BH formed at radius r, parameterized as tsuiv{r) = 
taurv,o{r /roY , we ask what is the average survival time 
(taurv) accounting for radial wandering. Noting that the 
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proper average is over the destruction rate, rather than the 
destruction time scale directly, we have: 

1 



(t 

\''surv/ 



(22) 



where rwand is the minimum radius from the SMBH to which 
the MSP-BHs have time to wander before being destroyed. 



Assuming nMSP-BH(r) 

(^surv ) 



(since Mns < Mbh), 



(23) 



drr^r °^^tsurv,o(''Ao)' 



For S = ctBH — 1 /2, the scaling appropriate for the dominant 
binary evaporation process in a Keplerian potential, and for 
our fiducial choice of qbh ~ 1.75, this yields 

1.25 



(^surv) ^ 1.25 



ro 



In 



/ Rn 



(24) 



7-) / I I ''surv,0 

^ ^max / \ ^wand / 

for rwand < i?max To evaluate the In (7?max/rwand) factor, 



we again exploit the loss cone scaling in equation (211: it 
corresponds to the solution of the equation ti^ = (tsurv) with 
Imin = ^wand, i.e. the miuimum radius that can be reached 
in the expected lifetime of the binary. Approximating the 



denominator of equation (21 1 as constant yields 



In 



I'lnax 
^ ^wand 

from which we obtain: 



R„ 



/{Rn 



1 1/2 



urv ''max J 

2.3 Gyr 



(isurv) ~ [1.5 Gyr X isurv(-Rmax)] ^ ■ 



(25) 
(26) 



Pulsar shut off: Although it does not destroy the NS-BH 
binary, if the neutron star stops emitting in the energy 
band of interest (e.g., in the radio), it is no longer an 
attractive observational target. The spindown time scale 
can be estimated as tspin = P/2P, where P is the spin 
period of the pulsar and P is its derivative. For MSPs, 
this time scale is typically a few Gyr. Since they are likely 
to continue emitting even after this time has elapsed, we 
assume that MSPs effectively live forever. Ordinary pulsars, 
on the other hand, have radio lifetimes < 100 Myr (e.g., 
Faucher-Giguere & Kaspi 2006 1 which are too short for 



a significant fraction to survive long enough to undergo 
an exchange interaction and be observed in a BH binary 
created as considered here. 

Age of the Galaxy: Regardless of the expected lifetime 
of the MSP-BH binary, the creation process cannot have 
operated for longer than the age of the Galaxy, tcai. 
We therefore always cap the survival time at this value, 
fiducially adopting toai ~ 10 Gyr. 



2.4 MSP-WD Progenitors in the GC 

The creation rate and survival time considerations above 
require an estimate of the starting number of MSP-WDs in 



order to yield a prediction for the number of MSP-BHs po- 
tentially observable today. As this is a complex population 
synthesis problem in itself, we do not attempt a detailed 
calculation, but instead base our results on simple physical 
scalings. The starting point for our estimate are the observed 
MSPs in dense globular clusters, such as Terzan 5 (Ter5), 
whose core mass density is similar to the density at the char- 
acteristic radius of the central BH cluster. Ter5 already has 
over 30 pulsars detected in the radio, most of them MSPs. 
Accounting for selection effects, including the fact that the 
cluster has only been probed down to a finite radio fiux, it 
could easily host 100 or more MSPs. Additionally, many of 
the globular cluster MSPs are found in binaries which fre- 



quently eclipse owing to ablation of the companion (e.g., Ta- 
vani||199T l, reducing the sensitivity of periodicity searches. 
The same MSPs would however no longer eclipse after ex- 
changing their companion for a stellar BH, and would there- 
fore become easier to detect. 

On theoretical grounds (e.g., Verbunt fc Hut|1987 l, the 
number of dynamically-formed binaries should scale with 
the collision number 



2 3-1 



(27) 



where po is the cluster density, is its characteristic radius, 
and CTo is its ID velocity dispersion. This expectation is 
well supported by observations showing that the collision 
numbers of globular clusters correlate well with the num- 



bers of close X-ray binaries and MSPs they host (Pooley 
|et al.|[2003l [Pooley fc Hut||2006| |Abdo et al.]|2010| ) In 
extrapolating equation (27 1 from globular clusters to the 



GC, other factors must however be taken into account: 

The number of NSs per unit mass: T he high kick ve- 
locities Wkick ~ 400 km s~^ of pulsars (e.g., Faucher-Giguere 
|fc Kaspi]|2006[ ) relative to the escape velocities Ve ~ 50 km 
s~ of massive globular clusters imply that only a small 
fraction of the NSs formed in globulars should be retained. 
Assuming that the birth kick velocity is Maxwellian in 
form, with (ukick) = 380 km s~^ (Faucher-Giguere & Kaspi 



2006 1, only 0.2% of the pulsars have Ukick < Ve = 50 km s~^. 
In practice, the retained fraction must be ~ 10% in order 
to explain the number of NSs observed in rich clusters, with 
the enhanced retention fraction most likely owing to the 
frequent occurrence of massive binary companions around 



the NS progenitors (e.g., Drukier 1996 Davies & Hansen 
11998'i Pfahl ct al.||2002[). Within the sphere of influence of 



the SMBH at the GC, the escape velocity is much larger. 



Ve{r) « 262 km s 



M 



4 X 106 Mq 



1/2 



0.5 pc 



-1/2 



(28) 

so that a fraction of order unity of the NSs born there could 
be retained, i.e. up to ~ 10 x as many per unit mass as in 
massive globulars. 

Recently, observations have furthermore suggested that 
the initial mass function (IMF) of stars at the GC may 
be top-heavy (Bartko et al. 20101. If that is the case, a 



larger fraction of the main sequence stars should produce 
NSs than in globulars. Let us assume, for illustration, 
that practically all stars with mass > 8 Mq produce a NS 
upon death, with a small fraction yielding a BH instead. 
Suppose further that the IMF covers stellar masses from 
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Mmin = 0.1 Mq to Mmax = 100 Mq . If the IMF has an 
unbroken Salpeter slope, dN/dm oc m' each 10'' Mq 
of star formation yields 7,400 NSs. If, on the other hand, 



the IMF is top-heavy with dN/dm 



throughout 



(as measured by Bartko et al. 2010 for massive stars), the 
same total mass of star formation produces 21,200 NSs, 
nearly 3x as many as the fiducial Salpeter case. In the 
intermediate scenario in which the IMF has a Salpeter slope 
below 8 Mq and -0.45 above, 16,300 NSs are produced, an 
enhancement of a factor of 2.2 over the pure Salpeter IMF. 
Combining a possibly more efficient NS production rate 
and a higher retention fraction, the GC could therefore be 
richer in NSs than globular clusters, per unit mass, by as 
much as a factor of ~ 20 — 30. 

Binary fraction: The dynamical interactions consid- 
ered in this work can only operate in the GC if binaries are 
present. There are currently almost no direct constraints 
on the binary fraction in the GC, precluding estimates at 
better than the order-of-magnitude level. Observations of 
hypervelocity stars (HVSs) (which in the classical model 
arise from the dynamical interaction of binaries with the 
central SMBH; Hills 1988), of the S-stars around Sgr A* 
(which may be the counterparts of HVSs captured by the 
SMBH; Gould & Quillen 2003, Ginsburg & Loeb 2006), 
and the detection of several X-ray binaries with projected 
separation <1 pc from Sgr A* (including one at < 0.1 pc; 
Muno et al. 2005a, b. Bower et al. 2005) however make a 
strong case that binaries do exist in non-negligible numbers 
in the GC. Assuming that ~ 10% of stars are members 
of a binary with semi-major axis ah ^ 0.3 AU, as in 
the Galactic disk (e.g., Duquennoy & Mayor 19911 and 



similar to what is inferred for globular clusters (e.g., Hut 



et al. 



1992 Rubenstein & Bailyn 19971, Yu & Tremaine 



(20031 predicted a HVS production rate roughly consistent 
with what is observed ( [Brown et al.|[2007[ ). It is therefore 
reasonable to assume that the GC is not severely depleted 
in binaries. 

Let /ns be the factor by which the GC retains more 
NSs per unit stellar mass than massive globulars, and let 
/bin be the ratio of the binary fractions at the GC and in 
globulars. We account, roughly, for the fact that binaries 
may be ionized down to smaller semi-major axes in the 
GC relative to globular clusters (as a result of the higher 
velocity dispersion) by adopting /bin = 0.5. For a binary 
semi-major axis distribution that is flat in the logarithm, 
the effect of truncating the upper end should in fact be 
modest on the total number of binaries. From the above 
considerations, we have: 



JVGC 



10 



MSP 



/ /ns 
V 20 



,/bin 
05" 



(29) 

where 'gl' stands for 'globular cluster' and 'GC stands for 
'Galactic center' as before. Here, A^msf s-nd A^'^'gp are the 
numbers of MSPs, assuming the same total stellar mass. 

For this comparison, we adopt the globular cluster pa- 
rameters ofTter5^p^ l^^jvl^pc^f^ 0.1 pc, af' = 12 
km s"^ ( [ Webbink| 1985) [Cohn et al.|2002[ ) . For the Galactic 
center, we assume a characteristic radius r^'~^ = 0.5 pc for 
the cluster of stellar BHs, and evaluate the stellar density 



and velocity dispersion at this radius. From equations ([T]) 



and 

f^o = 
yields: 



and for the GC model described in [2.1 we find 
139 km s"^ and = 3.3 x 10^ Mq pc'^ This 



■GC 
MSP 



12 



/ns 
20 



/bin 
0.5 



(30) 



Since Ter 5 and the central parsec both contain a stellar 
mass « 10® Mq , the above estimate implies that the central 
parsec could contain as many as 12 x more MSPs than Ter 
5, i.e., up to 1,200 ff Ter 5 hosts 100 of them. 

There are important uncertainties in extrapolating from 
globular clusters to the Galactic center. We therefore cau- 
tion the reader against over-interpreting this crude calcula- 
tion. Nevertheless, assuming that the MSP formation chan- 
nels operating in globulars operate in the GC as well (as 
supported by the observed overabundance of X-ray binaries 
in the central parsec; e.g., Muno et al. 20051, several hun- 



dreds to more than a thousand MSPs could populate the 
nuclear BH cluster today. We will adopt a nominal value 
A^MSP = 500 in the estimates that follow, approximately in 
the middle of the expected range. This number consistent 
with the upper limit of 90 ordinary pulsars in the central 



parsec reported by Macquart et al. (20101, as their observa- 



tions were insensitive to MSPs due to scattering pulse broad- 
ening at their observing frequency (15 GHz). Furthermore, 
globular clusters typically contain many more MSPs than or- 



dinary pulsars. As recently shown by Abazajian (20101, the 
7— ray fiux and spectrum measured by the Fermi Gamma- 
ray Space Telescop^ (i^ermi) toward the GC (over a larger, 
175 pc radius) is also consistent with a enhancement of 
MSPs in this region. 

For our fiducial parameters, we finally find 



J^MSP- 



f Nmsp 
V 500 



(31) 



2.5 Sensitivity to Model Parameters 

The calculations presented above assumed fiducial parame- 
ters for definiteness. In order to estimate the sensitivity of 
our results to the model parameters, we have repeated the 
calculations varying key ones, one at a time. Specifically, we 
have explored the variations listed in Table [l] in which we 
have varied each parameter by amounts representative of the 
expected range (for the case q:bh 7^ 1.75, the survival time 
is evaluated using a generalization of eq. ( [26| l). In addition 
to actual physical parameters, we have varied the Galaxy 
"age" icai as a test of the steady-state approximation. In all 
cases, the predicted number of surviving MSP-BHs is within 
a factor of 2 of the fiducial result in equation (311. At the 



order-of-magnitude level of the present work, our estimates 
are therefore robust to the details of the assumptions. 

The white dwarf mass, Mwd, warrants further discus- 
sion, as it also affects the predicted properties of the MSP- 
BH systems. In fact, the dimensionless ratio Mwd/Mbh was 
shown in §2.3| to determine both the factor by which the 
semi-major axis is multiplied during the exchange interac- 
tion and the eccentricity of the resulting MSP-BH. While 
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Table 1. Dependence of the Number of MSP-BHs Surviving in the Galactic Center Today on Model Parameters 



Parameter Varied 


^ "-min 


^max 








A'f-a-u 


^ ' BH 


'^bjHiin 


^b,max 


''Gal 


MSP — BHv^^^/ MSP^ 




pc 


pc 








A fr-, 
IVIQ 




AU 


AU 


dvr 




Fiducial 


0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


10 


6 


^min 


0.01 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


10 


7 


^min 


0.2 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


10 


6 


^max 


0.1 


1 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


10 


6 


"BH 


0.1 


0.5 


1.5 


1.3 


1 


10 


25,000 


0.01 


0.3 


10 


12 


CUMSP 


0.1 


0.5 


1.75 


1.5 


1 


10 


25,000 


0.01 


0.3 


10 


7 


MwD 


0.1 


0.5 


1.75 


1.3 


0.5 


10 


25,000 


0.01 


0.3 


10 


4 


MwD 


0.1 


0.5 


1.75 


1.3 


0.2 


10 


25,000 


0.01 


0.3 


10 


3 


Mbh 


0.1 


0.5 


1.75 


1.3 


1 


7 


25,000 


0.01 


0.3 


10 


5 


Mbh 


0.1 


0.5 


1.75 


1.3 


1 


13 


25,000 


0.01 


0.3 


10 


8 


A^bh 


0.1 


0.5 


1.75 


1.3 


1 


10 


15,000 


0.01 


0.3 


10 


5 


Nbh 


0.1 


0.5 


1.75 


1.3 


1 


10 


20,000 


0.01 


0.3 


10 


6 


^b,min 


0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.05 


0.3 


10 


6 


^b,min 


0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.1 


0.3 


10 


7 




0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


1 


10 


5 




0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


3 


10 


5 


*Gal 


0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


5 


6 


*Gal 


0.1 


0.5 


1.75 


1.3 


1 


10 


25,000 


0.01 


0.3 


2 


4 



Note. — All symbols are defined in the text. The fiducial parameter values are given in the top row and the model parameters are 
varied one at a time relative to the fiducial set, with the variation indicated in bold. All other parameters, such as the mass of the central 
massive black hole and the total mass distribution in the nuclear star cluster, are held fixed. 



the black hole mass also appears in this ratio, observations 
suggest that they have a narrower mass distribution (e.g., 
Ozel et al.|2010[). The fiducial choice of Mwd = 1 Mq, on the 



other end, is at the upper end of the observed white dwarf 
mass distribution, with values as low as Mwd ~ 0.1 M0 also 
occurring (even lower mass WD companions have been ob- 
served [e.g. Stairs 2004], but they are rare). Somewhat higher 
eccentricities than for the fiducial case are therefore possi- 
ble, e.g. e' ~ 0.98 for A/wd = 0.1 Mq, although the scatter 
in eccentricities resulting from exchange interactions is suffi- 
ciently large that this value would actually be perfectly com- 
patible with AfwD = 1 Mq as well (e.g., Sigurdsson & Phin- 



ney 



19931. In principle, MSP-BH semi-major axes larger 
than the ~ 3 AU expected from a Mbu/Mwo = 10 exchange 
with a binary with original semi-major axis flb.max ~ 0.3 AU 
are also possible for Mbh/ Mwd ^ 10, but these should be 
quickly evaporated (Fig. [T|. 



3 DISCUSSION 

3.1 Observational Prospects 

Our results suggest that our Galactic center should harbor 
several MSP-BH binaries, providing a strong motivations 
for focused searches in this direction. As no pulsar has so 
far been detected in the central parsec, it is important to 
consider whether finding MSP-BHs there is feasible. 

The time scale by which a radio pulse is broadened ow- 
ing to interstellar scattering is 



Tscat = 116 ms 



100 pc 



(10 GHz) 



(32) 



where Dscat is the distance of the eflective scattering screen 
from the GC. Combin ing all known tracers of ionized gas, 

133 



( 1998 1 found 73scat 



^gg" pc; we assume 



Lazio & Cordes 
i'scat ~ 100 pc. Sensitivity to a pulsed signal of period P 
requires Tscat ^ P- For a period P = 5 ms at the GC, the 
minimum observing frequency is therefore 22 GHz, while 33 
GHz is sufficient to detect a 1 ms source. Let us consider a 
hypothetical 30 GHz search, which would suffice to detect 
most MSPs if sufficiently deep. The minimum flux density 
for a pulsar detection is 



S„ 



She 



l^sys^iT'ccc. ~t~ ^sky) 



We 



(33) 



where Jbcam is a factor accounting for the reduction in sen- 
sitivity to pulsars located away from the center of the tele- 
scope beam, Tree is the receiver temperature on a cold sky, 
Tsky is the sky background temperature, G is the antenna 
gain, Np is the number of polarizations summed, Av is the 
receiver bandwidth, tint is the integration time, P is the 
pulse period. We is the effective pulse width, a is the signal- 
to-noise (S/N) detection threshold, and /3 is a constant ac- 
counting for various system losses (e.g., Dewey et al.|198"4 |. 
[Macquart et al.] ( |2010[ ) calculated the S/N expected as a 
function of observing frequency for a 10-hour observations 
with the GBT with a bandwidth Au = 800 MHz, for dif- 
ferent pulsar periods and spectral properties. We can use 



equation ( 33 1 to scale their results and explore the potential 



capabilities of future searches. 

[Macquart et al.| ( |2010[ ) predict a S/N = 1 for a 10-hour, 
30 GHz GBT integration with Af = 800 MHz, for a 5-ms 
pulsar with a typical spectral index a — —1.7 (5^ (x v") 
and flux density of 1 mjy at 1 GHz (corresponding to 
iiooo = 5ioood^ = 72.25 mJy kpc^ at the GC, d = 8.5 
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kpc), and S/N = 10 if the spectral index is instead -1. For 
reference, the observed spectral index distribution is well 
modeled by a normal distribution with mean —1.7 and dis- 
persion 0.35 (e.g., |Smits eraL][2009| . To get a sense of the 
limits of a dedicated GBT campaign with upgraded instru- 
mentation, we consider a 100-hr integration with a Av — 8 
GHz bandwidth insteacj^ in which case the same S/Ns are 
achieved for a pulsar with flux density lower by a factor of 
10. A futuristic Square Kilometer Arrajj^ (SKA) telescope 
would have an effective area larger by a factor of ~ 100, 
and hence would be sensitive to pulsars ~ 100 x fainter still. 
Thus, an upgraded GBT campaign could in principle detect 
a I/iooo ~ 7 mjy kpc^, a = — 1, 5-ms MSP at the GC with 
S/N — 10, while the SKA could in principle probe as faint 
as Z/iooo ~ 0.7 mJy kpc^ under the same conditions. For a 
more typical spectral index a = 1.7, the same S/N is reached 
for pulsars more luminous at 1 GHz by a factor of 10. For 
comparison, the least luminous MSPs currently detected 
in deep globular cluster observations have Liooo ~ 1 mJy 
kpc^ and the luminosity fun ction extends beyon d at least 
iiooo ~ 100 mJy kpc^ (e.g., [Hessels et al.|[2007[). A dedi- 



cated GBT search should therefore be capable of detecting 
MSPs at the GC, while the SKA could potentially probe a 
substantial fraction of the population. 

Fermi has demonstrated that MSPs are 7— ray sources 
(e.g., |Abdo et al.|2009 , Faucher-Giguere fc Loeb|2010| ). The 



7— ray energy band has the advantage that, unlike the ra- 
dio, it is not affected by interstellar scattering. However, 
even Fermi can only detect the very brightest, tip-of-the- 
iceberg 7— ray pulsars at the GC, so it is unlikely that a 
significant number of individual 7— ray MSPs will detected 
in the central parsec for the foreseeable future. 



3.2 Signatures of the Proposed Formation 
Channel 

Our calculations make definite predictions for the properties 
of the PSR-BHs formed via the proposed channel: 

(i) The system would be found within ~ 1 pc of Sgr A*, 
where the cluster of stellar BHs is assumed to be located. 

(ii) The pulsar would be recycled, with a period from ~ 1 
to a few tens of milliseconds, in order to have an original WD 
companion and to have a sufficiently long spin down time to 
remain emitting in the radio to the present day. Accordingly, 
the pulsar should have a low magnetic field, B < 10^" G. 

(iii) The MSP-BH binary would be relatively wide, as the 
semi-major axis is multiplied by a factor ~ Mbh/A/wd ^ 10 
during the exchange interaction. For an original MSP-WD 
semi-major axis distribution ranging from ~ 0.01 AU to ~ 
0.3 AU, the MSP-BHs should have semi-major axes ranging 
from ~ 0.1 to > 3 AU. 

(iv) The MSP-WD would be highly eccentric, 1 - e' ~ 
(MwD/MBH)(MBH/AfNs)''" (J 2.3 1. For our fiducial masses, 
e' ~ 0.8, but the eccentricity distribution will be broad. 



^ A spectrometer capable of such high-frequency, large band- 
width observations is already being developed for the GBT 



( http://www.gb.nrao.edu/gbsapp/ ). 
e- 



These predictions are based on the properties of the 
MSP-BHs immediately after their formation in exchange in- 
teractions. In principle, their internal properties could evolve 
by the time they are detected. Simple estimates however 
suggest that they should not change qualitatively. As dis- 
cussed in §2.3| the spin down time scales of recycled pul- 
sars are typically of order of a few Gyr. In the s oft r egime, 
the binary evaporation time tcvap oc a'^^ (eq. (16l), and 



so the binaries spend most of their lifetime with a semi- 
major axis of the order of its value at formation. In the hard 
regime, the binary semi-major axis is progressively reduced 
by hardening, but the time scale for this process is compa- 



rable to the exchange interaction time scale (e.g., Heggie & 
Hut|1993| l, which for MSP-BHs at the GC is > 10 Gyr (Fig 



Ij). These simple considerations are in good agreement with 
the Fokker-Planck modeling of binaries at the GC by |Hop-| 
( 2009 1, who found little evolution of the internal binary 



properties. Thus, "Heggie's law" (according to which soft bi- 
naries become softer and hard binaries become harder) holds 
but the time scales are too long for it to have a large impact. 

We use the cross section derived by Heggie fc Rasio] 
( 1996 ) to estimate the time it takes for encounters with other 
stars to change the eccentricity by an amount > Seo, start- 
ing with a value e. Averaging over a Maxwellian velocity 
distribution, 



te{> Seo) ~ 0.29 



0"3 



GnsMi23a'K 



M12M1 



1/3 



5e, 



2/3 



(34) 



e2/3(i_e2)i/3' 
where M12 = mi+m2 and M123 = mi+m2+m-t. Comparing 



with the ionization time scale given by equation ( 15 1, we find 



te(> Seo) 



.4/3 



9.69 



Se, 



2/3 



1/3- 



(35) 



http: / / www.skatelescope.org 



For our fiducial MSP-BH mass choices, e = 0.8 and Seo = 
0.1, te{> 5eo)/tic,n ~ 1.9, indicating that the eccentricity 
changes should also be modest over the lifetimes of the bina- 
ries. It is also noteworthy that even if eccentricity perturba- 
tions occurred on a shorter time scale, the systems would be 
asymptotically driven to a thermal distribution, /th(e) = 2e 
with mean {e)th ~ 0.67, so that the eccentricities would re- 
main large on average. 



3.3 Applications 

The discovery of a MSP-BH system at the GC would have 
profound implications. At this time, it would be the first 
known PSR-BH binary and would therefore critically in- 
form our theories of stellar evolution, especially in dense 
stellar environments. If it has the properties outlined above, 
it would convincingly originate in an exchange scenario sim- 
ilar to the one we have explored. It would therefore provide 
strong evidence for the existence of the predicted cluster of 
stellar BHs around Sgr A*, for which there is currently no 
observational confirmation. This would in turn have impor- 
tant implications for a number of other phenomena, includ- 
ing HVSs (e.g., [O^Leary fc Loebp OOS), gravitational wave 
sources (particularly for the extreme mass ratio events to 
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be detected by L/g^Dbut a lso for X/GCQ e.g., Barack fc 
Cutler] [20071 [O'Leary et"aL][2009l |Amar(>Seoane fc Preto 



20101, microlensing near Sgr A* (e.g., Chaname et al.|2001 



Alexander fc Loeb|2001 1 , and the orbital capture of stars by 



Sgr A* (e.g., [Alexander fc Livio|2004 | 



It would also for the first time provide an accurate clock 
orbiting a BH and therefore offer a unique probe of the 
spacetime in a BH potential (e.g.. Stairs 2003 Kramer 2007 1. 



While the MSP-BHs formed in this scenario would not be 
as tightly bound as the relativistic double NS systems from 
which the best constraints on gravity theories are currently 
derived (e.g., Kramer et al.|2006 Weisberg et al.|2010 1, cer- 
tain relativistic effects including the periastron shift and the 
Shapiro time delay can be measured even in relatively loose 
systems (e.g., [Taylor fc Dewey[p88| [Ryba fc TayIor][1991 



Narayan et al.|1991 1 . Measuring these effects would provide 
much more accurate measurements of the masses of stellar 
BHs than currently possible in X-ray binaries, in which as- 
sumptions regarding the inclination of the binary must be 
made. Another interesting possibility for nearly edge-on sys- 
tems would be to measure the gravitational lensing effects of 
the BH using pulse delays, from which the orientation and 



spin of the BH could be inferred (e.g., Narayan et al. 
Broderick fc Loeb|2006[ ). 



1991 



3.4 Caveats 

While our calculations provide straightforward predictions 
for where and how to find MSP-BH binaries, of their ex- 
pected properties, and suggest that they should be obser- 
vationally accessible, they are limited in some respects. As 
a first investigation of the proposed formation channel, we 
have focused on simple analytic estimates of the relevant 
physical processes, and modeled the GC as a steady-state 
background in which MSP-BHs are created and destroyed. 
In reality, the GC is a dynamical system in which compli- 
cated stellar evolution processes are ongoing. It would there- 
fore be desirable, as the capabilities become available in the 
future, to perform more detailed dynamical simulations that 
include the critical binary and stellar evolution effects. 

Some of the assumptions we have made have also not 
yet been empirically confirmed. In particular, our results rely 
on the presence of a central cluster of ~ 25, 000 stellar BHs 
induced by mass segregation, but this cluster may not be 
present if relaxation is too inefficient or if a binary massive 
black hole recently destroyed it (e.g., Merritt[|2010 1. While 
the disk morphology of the Milky Way argues against a ma- 
jor merger in the last ~ 10 Gyr, intermediate mass black 
holes brought in by dwarf galaxies or globular clusters can- 
nibalized by the bulge could also eject stellar BHs from the 
central cluster. Resonant relaxation, which we have mostly 
neglected, could furthermore be more efficient at depleting 
the stellar BH cluster or binaries near Sgr A* than we have 
assumed. It will be important to consider these issues in 
more detail in the future, although detecting MSP-BHs in 
the central parsec would actually inform these open ques- 
tions. 



http://lisa.nasa.gov 
* http://www.ligo.caltech.edu 



3.5 Conclusion 

We have shown that if dynamical processes analogous to 
those operating in dense globular clusters occur in the cen- 
tral parsec of our Milky Way, and if this region hosts a clus- 
ter of ~ 25, 000 stellar-mass BHs as predicted by mass segre- 
gation arguments, then MSP-BH binaries should be formed 
there in exchange interactions. Taking into account the much 
higher retention fraction of neutron stars in the Galactic 
center relative to globular clusters, as a result of the deeper 
potential well, several of these systems should survive to the 
present day. Simple scaling arguments indicate that some 
of these MSP-BHs might be detectable with existing radio 
telescopes, and that a substantial fraction of the population 
should be accessible to the Square Kilometer Array. Our 
predictions therefore suggest a definite roadmap to the de- 
tection of the first pulsar-black hole binary, by singling out 
a small region of the Galaxy where many might reside. In 
light of the remarkable potential payoffs of such a discovery, 
it is therefore clear that focused observational searches and 
further theoretical studies are warranted. 
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